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Estimate of the charming penguin contributions to B 
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Abstract. We consider the problem of factorization in B decays and illustrate the calculation of non- 
factorizable contributions employing the QCD light-cone sum rule method. We present a more detailed 
calculation of the "charming penguin" contributions as a potential source of the substantial nonfactorizable 
0(l/rri),) effects in the B — > nn decay. Although the predicted corrections are not sizable, by calculating 
the CP asymmetry we illustrate how such corrections can accumulate to a visible effect. In conclusion, 
nonfactorizable contributions in nonleptonic B decays into charmonium are briefly discussed. 
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1 Nonleptonic B decays and factorization 

For a long time the naive factorization method, in which 
the matrix element of the four-quark operator is approxi- 
mated by the product of two matrix elements of the bilin- 
ear quark currents, was considered as a sufficiently precise 
tool for estimating matrix elements emerging in the am- 
plitude of nonleptonic weak B decays. Nowadays, in order 
to make real use of the already very precise experimental 
data, we are forced to provide a more accurate estimate 
of nonleptonic decays, in particular of the nonperturba- 
tive part of the decay amplitude. Therefore, the question 
about the applicability of the factorization and the size of 
nonfactorizable corrections naturally emerged. The ques- 
tion was particularly raised in the work [Tj, where it was 
argued that in the charmless B decays there could exist 
large 0(Aqcd /nib) corrections and large strong phases 
coming from the " charming penguins" . 

There are several models which one can apply for the 
calculation of matrix elements of B-meson weak decays 
beyond the naive factorization . By using the QCD 

factorization approach [3] one can show that the exclusive 
_B-decay amplitude in the mj, — > oo limit can be expressed 
in terms of the factorizable part and the calculable 0{a s ) 
nonfactorizable correction. However, because of the argu- 
ments given above, the nonfactorizable subleading effects 
in the decay amplitude, suppressed by inverse powers of 
mb, could be important and have to be investigated. Es- 
timates of nonfactorizable contributions in B decays, in- 
cluding the power-suppressed 0(l/m,b) contributions, can 
be obtained fl] using the method of QCD light-cone sum 
rules (LCSR). In particular, nonfactorizable contributions 
to nonleptonic decays such as B — > rnr and B — > ttK 
become interesting for a more precise constraint on the 
7 = arg(V u b) angle of the CKM matrix. In connection 
with this problem, it is worth mentioning that there are 



also several strategies used to determine the 7 angle from 
B — ► 7T7r and B — > ttK based on the isospin and SU(3) 
relations. Unfortunately, the theoretical accuracy of these 
relations is limited and it has to be improved by calculat- 
ing the SU(3) breaking effects, which can also be addressed 
by the LCSR method 0. 



2 Matrix elements for B 



TTTY 



from LCSR 



The LCSR expression for the B — + 7T7t hadronic matrix 
element of the Oi operator of the weak Hamiltonian is 
derived by the procedure presented in detail in fJ]|H] • One 
starts by introducing the correlation function 

^\T{3 { :kv)om3f\m^-^)) 

= {p-k) a F^\ Sl ,s 2l P 2 ) + ..., (1) 

where j£g = ujaj^d and = imjrysd are the quark 
currents interpolating the pion and the B meson, respec- 
tively. By employing the dispersion relation technique and 
by assuming the quark-hadron duality, we can write the 
LCSR expression for the hadronic matrix element A^ 0i ^ (_B° 
7T+7T-) = (Tr-(p)ir+(-q)\Oi\B d {p - q)) as 0] 
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where M\ and M 2 are the Borel parameters in the pion 
and -B-meson channels, respectively. The parameter Sq 
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(sq) is the effective threshold parameter of the perturba- 
tive continuum in the pion (B-meson) channel. In the sum 
rule J3J the finite corrections are taken into account, 
but numerically very small corrections of order sj /rn? B are 
neglected. 
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(a) (b) 

Fig. 2. Contributions in the penguin topology: (a) an exam- 
ple of the 0(a a ) nonfactorizable penguin amplitude; (b) an 
example of the chirally enhanced twist-3 contribution to the 
penguin amplitude. The square stands for ' c and Oi-n oper- 
ators. The leading contribution from the Og, operator proceeds 
without the quark loop |8I7| . 



(a) 



(b) 




(c) 

Fig. 1. Contributions in the emission topology: (a) factoriz- 
able contribution from the 0" operator; (b) an example of the 
0(ct s ) correction; (c) an example of the nonfactorizable dia- 
gram from the (D\ operator. 



soft (low-virtuality) gluons coupled to the c-quark loop 
□ in the sum rule approach manifest as multiparticle 
DA's and are therefore suppressed at least by 0(a s /mt). 
Therefore, the nonfactorizable 0(l/m b ) corrections from 
penguin loops are mainly of perturbative origin, and both 
contributions from Fig. 2. generate the strong rescattering 
phases in B — > tttt perturbatively by the well-known BSS 
mechanism gjj. 
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3 CP asymmetry in B 
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2.1 Emission topology 

The leading contributions in the emission topology are 
shown in Fig.l. The factorizable part (Fig. la) stems from 
the contribution of the leading operator 0" = (dr^u)(ur^b) 
mainly reproducing the naive factorization result in terms 
of the light-cone calculated f v and Fb^tt form factor. The 
hard corrections (Fig. lb) have not yet been addressed in 
LCSR, but the corresponding contributions calculated in 
the m b — -> oo limit in QCD factorization .'5 appeared 
to be small. The soft nonfactorizable contribution due to 
the Of operator (which is the color-octet part of = 
1 /N c Oy + 2<D\ and expresses the exchange of a soft gluon 
between two pions (Fig.lc)) was shown 0] to be small, al- 
though of the same size as the 0(a s ) correction mentioned 
above. 



2.2 Penguin topology 

Types of the dominant penguin diagrams are shown in 
Fig. 2 [7|. The main effect which we calculate arises from 
the c-quark loop annihilation into a hard gluon (Fig. 2a). 
In addition, there is the quark-condensate contribution 
(Fig. 2b) which after a detailed analysis appeared to be a 
natural upper limit to all neglected contributions of multi- 
particle distribution amplitudes (DA's). The effects of the 



Penguin contributions appear to produce a notable effect 
in the direct CP asymmetry and we take the CP asym- 
metry as a testing ground for the influence of the l/m& 
corrections in the charming penguin contributions. Fol- 
lowing jjj , we concentrate on the direct CP asymmetry in 
t + 7t~ decay, which is given as 



the B° d 



-,dir 

X CP 



= (i-|£l + 



(3) 



where £ = e -M0+~t)(i + R e ^)/(1 + Re^) and R = 
—P/(RbT). Here T is the contribution to the B — > tttt 
amplitude proportional to V u bV* d = \V u bV* d \e~ n . It con- 
tains the tree amplitude, the penguin-loop contractions of 
the current-current operators 0" 2 > an d also the VubV* d 
proportional penguin O3-6 operator contractions. The re- 
maining contributions, being proportional to V c bV* d , are 
included in P. The penguin- loop contractions of the current- 
current operators 0\ 2 represent the main contribution to 
this part. The factor R b = \V ub \\V ud \ / (\V cb \\V cd \) is the 
ratio of the CKM matrix elements. 

Both T and P amplitudes have strong phases; there- 
fore, we have T = \T\e l5T and P = \P\e l5p and the CP 
asymmetry for B d — > 7r + 7r can be written as 



— 2\R\ sin(<5p — St) sinj 
l-2\R\ cos(5 P -5 T ) cos 7 +|i?| 2 



(4) 



All contributions shown in Figs.l and 2 are calculated in 
LCSR at finite m&. We also include the LCSR result for the 
gluonic penguin contribution of the dipole operator Og, g 
[5] . The electroweak penguin contributions to B — > tttt are 
color-suppressed and negligible. 
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Fig. 3. Direct CP asymmetry in 73° — » tt + iv~ as a function of 
the CKM angle 7. The uppermost curve is the result obtained 
for m,b — > 00. The dark region is the LCSR result, with all 
uncertainties from the method included (uncertainties in the 
CKM matrix elements are not taken into account). The light 
region shows the deviation from the mi, — > 00 limit result. 

The hard 0(a s ) corrections to T and P amplitudes are 
known in the mj — > 00 limit from QCD factorization 3 . 
We have examined the influence of these contributions to 
the phases St and Sp.lt appears that they are highly sup- 
pressed in comparison with the phases emerging from the 
penguin-loop contributions. Therefore, we have neglected 
0(a s ) corrections in J3J. 

In Fig. we show function of 7, calculated 

by using the penguin contributions estimated from LCSR 
at finite nib (dark region) and compare the result in the 
infinite-mass limit that agrees with the QCD factoriza- 
tion prediction [3] (the uppermost curve). Both results 
are taken at the same scale /i;, ~ rrib/2 as used in LCSR. 

The prediction shown in Fig. 3 is not final, annihilation 
effects are missing and the uncertainty in the CKM matrix 
elements is not taken into account either. However, the 
figure nicely illustrates the size of 0(l/rrib) corrections and 
the difference between the results obtained at the finite mb 
and in the nib — > 00 limit (the light region in Fig. 3). 

4 Conclusion 

We have discussed the factorization in nonleptonic B me- 
son decays and have presented the LCSR calculation of the 
B — > 7T7T decay. It has been shown that the leading contri- 
bution factorizes, while the corrections beyond the factor- 
ization can be systematically approached. In the m — > 00 
limit, our result agrees with the QCD factorization pre- 
diction, while at finite nib we have found 0(a s /rrib) effects 
which are numerically small, but accumulate to a sizable 
correction in the direct CP asymmetry. Large charming- 
penguin contributions per se are not predicted by this 
model. 

Numerically, the nonfactorizable corrections in B — > 
7T7T are not large, but the situation in the nonleptonic 
color-suppressed decays seems to be somewhat different. 



Recent measurements of the color-suppressed B — > D^ tt' 
decays, B — > J/ipK and particularly B — > XcjK decays 
show large discrepancy with the naive factorization predic- 
tion and provide clear evidence for large nonfactorizable 
contributions. The B — > XcoK and B —> Xc2K decays are 
particularly interesting because their branching ratios pre- 
dicted in the naive factorization are exactly zero, and the 
measurements yield branching ratios ~ O(10~ 4 ) which are 
therefore comparable with BR(B — > J/ipK). This should 
not come out as a surprise, because the LCSR calculation 
shows the existence of large nonfactorizable corrections in 
B — > J/tfjK of the order of 70% The corresponding 
large corrections are also expected for B — > XcjK decays. 
Unfortunately, these large corrections still appear to be in- 
sufficient to reproduce the data, because the discrepancy 
between the the theory and the experiment is much larger 
(in the case of B — » J/tpK, the data and the LCSR im- 
proved prediction still differ by a factor of two) and at the 
moment it is not clear how these discrepancies could be 
reduced. 
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